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Abstract: In this paper we shall analyse the structure of the Cauchy Problem (CP briefly) for General 
Relativity (GR briefly) by applying the theory of first order symmetric hyperbolic systems. The role of 
harmonic coordinates will be discussed. 


1. Introduction 

The Cauchy Problem (CP) for standard General Relativity (GR) has been studied in several 
papers, from numerical viewpoint (see, for example, [1], [2], [3]), as well as from an analytical 
viewpoint (see [4], [5], [6]). It is at the basis of all applications of numerical gravity as well as 
at the basis of the physical interpretation of GR theory. 

The analysis of the CP for GR is essentially based on Arnowitt, Deser, and Misner (ADM) 
seminal paper; see [ 7 ] and also [ 8 ], [ 9 ], [ 10 ]. The ADM decomposition is introduced and 
particular coordinate systems, called harmonic coordinates, are used to study the conditions 
under which one has a well-posed CP. 

In this paper we shall review such a procedure aiming to clarify the role of the choice of 
harmonic coordinates. In fact, in GR one is analysing covariant equations, the Cauchy theorem 
is stated in terms of the principal symbol of the differential operator which is defined intrinsically 
(as we shall show also in GR in which the equations are non-linear but quasi-linear and the 
principal symbol depends on fields) and hence it does not depend on the choice of coordinates. 
This is a first step in a project aiming to characterise integrability and constraints (as done in 
a Hamiltonian framework) though in a Lagrangian setting. 

The ADM decompositions transforms the original (vacuum) Einstein field equations (which are 
10 PDE in dim(M) = 4) in a new system formed by six equations, which will be a hyperbolic 
system, and four constraint equations. Only the first system of 6 equations is necessary to 
define a CP but we have to keep in mind that the second one is important as much as the first 
one for the physical problem and it constrains the allowed initial conditions. 

ADM decomposition also allows us to define a parameter r which represents the evolution pa¬ 
rameter of the system. In other words, from a mathematical point of view, ADM decomposition 
is a choice of a bundle structure {M, R, r, E) over the spacetime M, namely: 


( 1 . 1 ) 


r: M — >• R. 


The standard fiber E is a model for the isochronous space submanisolds Et^ := r ^(to) C M. 


The ADM transformations are the transformations which preserve the bundle structure: 



( 1 . 2 ) 


where are fibered coordinates over M (see [ 10 ], for further details). 
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On the spacetime M, one can restrict to Lorentzian metrics g for which the fibers T,t„ are 
space-like submanifolds. One can then decompose the metric tensor (or its inverse g^'^) in 
the following way: 


= 


f -N^ + |iV|2 

I Nj 



/ -iV-2 N-'^N^ \ 

\N-^N^ Y^-N-^N^Ni) 


(1.3) 


where TV is a spatial scalar field called the lapse, iV' a spatial vector called the shift and 7 y a 
3-Euclidean metric called the induced metric, defined on E, with respect to transformations on 
S (see [10]). 

Let ft be the (future directed) unit vector g-orthogonal to Ej and a basis for vectors tangent 
to the fibers Ej^. Then the lapse and shift are defined by the relation 


do = Nfi+N^ei 


(1.4) 


Ricci tensor, written in these new fields and in the frame (n, Ci), is: 

^00 = -^ {^oK - D,D^N) - 

< Roj = Di{Kj‘-6‘K)=R^, (1.5) 

^ ^ i7JlSoK^.^ - D.D,N) + ^R,, + KK,, 

where Jq = dg — £j;^, Kij is the extrinsic curvature (namely iVy = and is the 

covariant derivative with respect to the affine connection ^Th induced by 77 . 

The Ricci scalar is: 

i? = A (^SoK - DiD^N) +^R + K^ + ( 1 . 6 ) 


Einstein equations do not determine the evolution of N and N'^, i.e. they are not dynamical 
fields, so they can be arbitrarily chosen. We will set hereafter, for sake of the simplicity, TV = 1 


and W = 0 . 

With these choices we have that the metric tensor and its inverse become: 


9ri' = 


-1 0 

0 7y 




-1 0 

0 7^ 


(1.7) 


and the ADM evolution Einstein equations Rj^ = 0 become: 

ki{y}{,™}9o5oV’" + «0 (1.8) 

where we set 




(1.9) 


and: 


B 


kl _ kl 

{ij}{mn} ~ I 


'yi{m^n)j ^mn 


5^5^] 




■l^(l 




( 1 . 10 ) 


where the symbol « means modulo lower order terms with respect to the derivative degree. See 
[10] for further details. 
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Let us stress that the other Einstein equations i?oj = 0 and Jigg = 0 do not contain second order 
time derivatives (one can eliminate the term SgK by summing with the trace of the equation 
BLji = 0) and are hence to be interpreted as constraints on the allowed initial conditions. We 
shall not consider these constraints here (which contain information about the canonical analysis 
of the system; see [11]). 

Next section is devoted to state the Cauchy Problem for the evolution PDE (1.8). 


2. The Cauchy problem for a PDE 

First of all, let us remind that we are not analysing a generic PDE but quasi-linear systems 
which come from a Lagrangian. This means that we have a Lagrangian L = L{x^,'if^)dS 
and we obtain the Euler-Lagrangian equation by standard action variation. After the action 

Ad[(^]= f L{y{x),dy{x))dS (2.1) 

JD 

has been varied (with fixed boundary conditions) we obtain: 

SAdIu] = J aj{x,y)5y’dS = 0 (2.2) 

and the relative equations are 


oi.j{x,y)=Q. (2.3) 

It is clearly that if equations of motion come from a Lagrangian (by using the principle of 
least action) they will live in the dual space of the fields deformations . In general, by the 
geometrical framework for variational calculus (see e.g. [ 12 ]), one can show that Euler-Lagrange 
equations are described by a (vertical) bundle morphism 

E : j2Lor(M) y*(Lor(M)) (g) A^{M)) (2.4) 

which exactly expresses this remark. 

Let us briefly review the CP for a first order PDE. Although GR has second order equations let 
us introduce it for first order and then we shall extend the results to the second order system. 
Let us also stress that Einstein equation are quasi-linear, namely they are linear in the highest 
derivative terms. In this view, a first order PDE analogous to the one for standard GR for us 
is written as follows: 

oiij{x, y)dgy’ + a}j(x, y)diy^ + ji{x, y) = 0. (2.5) 

In the same way, in order to define a GP, we have to define an initial condition, namely: 

aij{x, y)dgy^ + a]j{x, y)diy^ + 7 /(a;, y) = 0 
y\0,x^)=f\x^) 

where f'^(x^) is called initial conditions or initial data. 
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We are ready to state the existence and uniqueness theorem: 

Theorem (2.7): Let ajj be a non-degenerate positive-deEnite bilinear form and let a\j be 
symmetric in indices IJ for all i, then under these hypotheses, given the initial data in H^(M), 
with k > ^ + 1, the existence and uniqueness is ensured in an open interval / C K and in a 
suitable Sobolev space, which depends on the regularity of the initial data (see [6] j. 


We stress that, under the above conditions, if we have smooth, we will have a smooth 
solution defined on / subset of the real line that contains t = 0. In other words, we imposed 
initial conditions and the equations determine a unique solution in a neighbourhood of the 
Cauchy surface So- 

Overviewing technical details, we can notice that two aspects are involved in the theorem 
above: an analytical and an algebraic condition. Although they are both important we shall 
focus on the algebraic one, also in view of the fact that in most contexts physical fields are 
chosen to be smooth. 

Let us remark that the theorem above states that the well-poseness of the CP is subjected to 
some algebraic conditions of the coefficients appearing in the differential operator. 

Now, we can consider the case of second order quasi-linear systems. As done in (2.5) we define 
a second order system as follows: 

aij{x, y)dooy^ - a}j{x, y)doiy^ - dlj{x, y)dijy^ + yi{x, y, dy) = 0. (2.8) 

and its associated CP (hereafter we drop the coefficients dependence): 
aijdmy’ - a\jdQ,y'’ - ajjd.jy’ « 0 . 

y\Q,x^) = f’{x^), ao/(0,x*) = /(x*). 

which will be called hereafter CP2. 


Our goal is to transform a second order PDE in a first order system, by introducing auxiliary 
fields. Inspired by the method used for ODE we can define the following new fields: 

{v'^ = , , 

K = (2.10) 

I < = djy^ 

However, if we wish to consider these equations as part of the original system one should notice 
that they do not correspond to differential operators with values in the dual space of field 
variations as it was for the original equation. Then one should introduce some suitable bilinear 
forms to write them equivalently in the form: 
i PiJ (doy^ — V^) = 0 

S ^ Pu K - djV ) = 0 (2-11) 

\ - d,v^) =0 ) 


for some invertible coefficients jdjj and /3)j. 

Let us remark that the equation vj = djy^ contains no time derivative and as such is a 
constraint on initial conditions and will not contribute to the CP. 


Then the equation (2.8) can be written in terms of the new fields {y^,v^,v(), so that the CP 
can be recast in the following form: 

(dijdoy ’« 0 

I aijdoV^ - a}jd,V^ - a)^jd,vj « 0 

yJiO,P) = f(P), v\Q,P)=y\P), vfiO,P)=d,f(P) 
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together with the constraint diy'’ = vf. We can write this system in the block-matrix form as 
follows: 


/pij 0 0 \ /y’X /O 0 

I 0 ajj 0 j i9o I 1 — I 0 a\j 

Vo 0 Vt-// Vo Pfj 


kj 

^ij 

0 



0 


(2.13) 


This is a first order CP so the theorem (2.7) applies to it. One has existence and uniqueness 
of solutions if the first matrix is symmetric, non-degenerate, positive-definite and the second is 
symmetric. 

We already know that ajj is non-degenerate and positive-definite. If also fjjj and /3)j are 
non-degenerate and positive-definite then the whole matrix is. Since we are free to choose /3/j 
as we wish (provided that the choice is non-degenerate and positive-dehnite) we can fix it as 
Pij = ajj, which is automatically a good choice. 

For the second matrix to be symmetric aJj must be symmetric and one must have (see Ap¬ 
pendix C) 

= (2.14) 


Thus the block Pfj (and as a consequence of this choice the coefficient Ojj) must be symmetric 
in (/J) and non-degenerate positive-definite. 

We can rewrite the original system as: 



0 




, / 

aij 

0 

1 1 


- 

0 



Kvj) 

' \ 


0 0 


U a 
0 a 


ij 


ij 


kj 

^ij 

0 


\ 


dk 


1 

Kvj / 


(2.15) 


V>{0,x^)=g-'ix^), vi{Q,x^)= hi{x^) 


which, together with the constraint diy^ = vf, is called the reduced CP or CPI for short. 

Let us remark that once again, also for second order operators, the well-poseness of the CP 
is subjected to algebraic requirements. Unlike for first order operators symmetry is no longer 
enough and one needs to require that the coefficient is also positive-definite. 

Now we have to show that the original CP2 (namely (2.9)) is dynamically equivalent to the 
reduced CP (2.15), namely CPI. 

Obviously, if we have a solution y'’{t,x^) of CP2 then: 


{y\v^ := doy\vf := 9,/) 


(2.16) 


is a solution of the reduced CPI. In fact, one immediately has that 

aij{v^ - doy-^) = 0 alj{dovj - djV’) = 0 (2.17) 

while the second order equation can be recast as 

aijdov^ — a\jdiV'^ — aj/jdivj « 0 (2-18) 

Thus the equations of CPI are satisfied, the constraint is satisfied (since we defined v/ := diy^) 
and the constraint, evaluated at t = 0, shows that hj(x) = djf’^{x) are the only initial conditions 
compatible with the constraint. Thus CPI holds true. 
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Viceversa, we need to prove that, given a solution {y\t,x),v'\t,x),vl{t,x)) of CPI which 
satisfies the constraint, then y\t,x) is also a solution of CP2. One has that = doy^, by 
first equation in CPI, and that dovj = djV^ = d^y^, by the third equations. Thus, one has 
vj = djU'^ + k{x) though, because of the constraint, k{x) = 0. Accordingly, one also has 

u/ = = 9o/ (2.19) 

with which the second equation of CPI implies the second order equation. 

One obviously has that y^{Q,x) = f'’{x). Now, the first equation of (2.15) tells us that: 

d()y\t, x) = v\t, a;) d^y’{Q, x) = u-^(0, x) = g’’{x) (2.20) 

which is the second initial condition and CP2 holds true. 

Obviously the correspondence 

{y\t,x)) ^ (/,z;^:=ao/,u/:=a,/) (2.21) 

sends solutions of CP2 into solutions of CPI and is a bijection, proving dynamical equivalence. 

We have eventually to prove that the if constraint is satisfied at the initial time, then it will 
be satisfied at all time. For, let us define the quantity 

fc/ := vj - (2.22) 


Since we know that doy' = v'’, then: 


dokj = dovj - dojy'^ = d^vj - djV'^ = 0 . 


(2.23) 


where the last equality holds true by the last equation of CPI. 


In particular, by imposing the constrain is t = 0, we obtain kj (0, x) = 0 and then we have the 
following Cauchy problem: 


otLjdokj = 0 

kj (0, a;) = 0 


(2.24) 


which has a solution kj (t, x*) = 0, which is unique since the system is symmetric hyperbolic. 
This means yj{t,x) = djy^{t,x) at any time and the constraint is satisfied at all times. 


Next section is devoted to apply the theory of PDFs developed above to the GR case. 


3. Harmonic coordinates 

We can start to reduce the system (1.8) of second order to a first order one. We have to define 
some new fields, as done in the previous section: 

( 3 . 1 , 

1 s;' := diY’ 

so that the system can be rewritten in the following way (as done in the previous section): 

f - AjyjpmjE'™ = 0 

I « 0 (3.2) 

I rikl a \Mnn rykl a \Mnn _ n 

I ^{ij}{mn}^0^k - ^ 
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with the constraint 


(3.3) 


sf = da^^. 


In matrix form the system above can be written as follows: 

'0 0 0 



0 



^ ym \ 

1 / 

0 


0 

1 ao 

1 

- 

0 

0 

TDkl i 



' \ 


0 

0 0 
0 


Qkl 

0 


dk 

f rylm \ 
Y^mn 1 

1 

i y-mn / 


0 . 


(3.4) 


We want to verify whether the hrst order system (3.4) is symmetric hyperbolic. We have two 
conditions that have to be satisfied, first one is that matrix 


B = 


0 
0 

' 0 }{mn} 


0 

Qkl 


(3.5) 


has to be symmetric which is automatically satisfied by construction. Second one is that matrix 


'A 


A = 


0 

0 


A 


0 


0 
0 

okl 


(3.6) 


has to be symmetric and positive-definite. One can easily show that the block yljyjpm} is 
positive-definite and symmetric. On the contrary symmetric with respect to the 

exchange of the pairs {ij}{lm}. 

In fact, it is easy to see that its antisymmetric part is: 


r)kl 




(3.7) 


and it does not generically vanish. Let us also remark that it is a tensor on the space manifold 
S so that one cannot hope it will vanish in any (spatial) coordinate system. 

Let us introduce a new coordinates system, called (spatial) harmonic coordinates, and defined 
by the following conditions: 

3^1 = ^ 0 ( 3 . 8 ) 

where are the Christoffel symbols of jij. It is easy to prove that harmonic coordinates 
always exist. The condition (3.8) is equivalent to: 


= 7“'9^7^fc 


which in turn implies 


Then we have 


zljndida^'^ « -iijdidni 


nj 


so that, in harmonic coordinates, the coefficient takes the form: 

-®{y}{mn} “ 7 'li{mln)j — 7 ^{ij}{mn} 


(3.9) 

(3.10) 

(3.11) 

(3.12) 


which is symmetric (as well as automatically non-degenerate positive-definite since ^{y}{mn} is- 
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In other words, in harmonic coordinates the antisymmetric part of the full operator becomes 
lower order, it contributes to lower order tail and consequently the operator becomes symmetric 
hyperbolic and CP is well-posed. 

Let us remark that there is no contradiction between B^{ij}{mn}] being a tensor and it becom¬ 
ing symmetric in harmonic coordinates since what we are really saying is that in harmonic 
coordinates one has 

^ 0 (3.13) 

not that the coefficients ^[pj}{mn}] become zero. In fact, the equations (3.13) are not covariant 
and they can be satisfied in particular coordinate systems (e.g. harmonic coordinates) without 
being satisfied in others. 

Now that we have seen that harmonic coordinates make evolution equation to be symmetric 
hyperbolic (though of course they spoil general covariance), we shall show that coordinates more 
general than harmonic coordinates exist for which the evolution equations are still symmetric 
hyperbolic. 

Since we have already found the antisymmetric part (3.7), we can directly impose that 

^ 0 (3.14) 

without imposing conditions on first derivatives as in (3.9). 

This condition is weaker than harmonic coordinates conditions. Obviously, harmonic coor¬ 
dinates imply (3.14). Viceversa, we shall see that there exists a coordinate system in which 
evolution equations become symmetric hyperbolic and this system is not harmonic; see [ 13 ]. 
In fact, the condition (3.14) is also satisfied if one simply has 

= 0 (3.15) 

Indeed, if the inverse metric takes a linear form in coordinates, e.g. 

/I 0 0 X 

7^^ = 0 1 0 (3.16) 

\0 0 1-ft/ 

then (3.15) vanishes, while P^ = 7b3r'y does not. 


4. Bianchi identities and constraints 

In this Section we shall see how one can use projectors (see appendix A) to prove that evolution 
preserves the constraints 

U {m - D,D^N) - = 0 

\ M, := Di {Kj‘ - 5\K) = 0 

due to (contracted) Bianchi identities. In fact, Bianchi identities are 


= 0 
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(4.2) 



which by projection provide the following conditions 


= 0 {5oH = -D,{NM^) - MW,N - 2NKH 

= 0 \ SoM^ = -D\NH) - 2NM^Kj,^ - NKM^ - HD^N 


By setting, as done above, = 1 and N = 0 and we obtain: 


doH « -S^dkM^ ^ ^ ] 

^ -S^dkH ^ VO 7m 7 °\M^) 


fo 6! 
\Sf 0 



(4.4) 


which is symmetric hyperbolic. Then, a unique solution exists for any initial condition. If we 
set H|t=o = 0 Slid ^\t=o ~ initial condition, then H = M* = 0 at alH G / is a solution 

and since equations are symmetric hyperbolic it is the only solution. Thus the constraints are 
preserved by Bianchi identities. 


5. Conclusions and perspectives 

We have seen that (vacuum) Einstein equations can be split into an (elliptic) system of con¬ 
straints (4 equations in dim(M) = 4) and an evolutionary system (6 equations in dim(M) = 4). 

The evolutionary part is not symmetric hyperbolic in general. However, one can split the 
evolutionary part into a symmetric hyperbolic equation and a further constraint, which is non¬ 
covariant with respect to change of coordinates on the spatial manifold E. Hence one can hnd 
spatial coordinates for which the antisymmetric part vanishes and, in that coordinates, solving 
the symmetric evolutionary part of equations. 

The evolutionary equation is covariant with respect to change of coordinates on S, thus if 
a unique solution is found in a coordinate system, then a solution is found in any coordinate 
system. This is not strange after all since being symmetric hyperbolic is a sufhcient (not 
necessary) condition for solving CP. 

Thus there exists a solution to the evolutionary part for any initial condition. However, not all 
initial conditions are the same: there are initial conditions which satisfy the elliptic constraints 
as well as initial conditions which do not. For any initial condition which satisfies the elliptic 
constraints one can find a spatial metric 7 y (t, x), which together with a choice of the lapse N 
and shift fields N (which in fact fixes the ADM foliation), defines a global Lorentzian metric g 
which solves original Einstein equations. 

This is more or less well known since it is the basis for numerical gravity (see [1], [2], [3]). 
However, a detailed analysis allows to clarify some of the details which are relatively less well 
known and to draw some conclusions which, to the best of our knowledge, are new. 

First we clarify the role of (spatial) harmonic coordinates. The coefficients A{ij}{mn} s-nd 
^{ij}{mn} ADM splitting of Einstein equations (i.e. their principal symbols) are spatial tensors. 
Hence, one can prove that changing the coordinates will not make them symmetric if they are 
not in the first place. 

This originally appeared as an issue to us: how can covariant CP be well posed in a coordinate 
system and not well posed in another? The solution is quite simple: the evolutionary system is 
well posed in any coordinate system, just in some coordinate system one can use the theorem 
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about symmetric hyperbolic systems while in other coordinates the same system is not symmet¬ 
ric hyperbolic (though the CP is well posed anyway). In other words, it is the condition of being 
symmetric hyperbolic that is not covariant. Nevertheless, the Cauchy theorem for symmetric 
hyperbolic systems, despite being not covariant, seems to be enough to deal with standard GR 
in full generality. 

As a first issue then, one should define an integrable evolutionary system as one which is 
symmetric hyperbolic in at least some coordinate system. 

Changing coordinates does not make the system symmetric hyperbolic. It simply makes their 
antisymmetric part identically satisfied. In other words, the antisymmetric part of the coefficient 
B^{ij}{mn}\ remains non-zero. However the equation from it 

^Ky}{mn}l^«7™ ~ 0 (5-1) 

becomes identically satisfied. This is possible precisely because this equation is non-covariant. 

Moreover, in view of a generalisation, we can stress that the fact that the original Einstein 
equations come from a variational principle play a fundamental role, in particular in GR when 
a metric on fields (namely, is unknown since it depends on the unknown field to 

be determined 7 y . In particular, one has as many equations as fields. If k equations are 
constraints then it is reasonable to expect that k fields will be left undetermined, eventually 
spoiling uniqueness as we know it must be in view of the hole argument (see [ 14 ], [ 15 ], [ 16 ]). 
Moreover, constraints are responsible for the fact that the system is overdetermined. This 
simply accounts for the fact that equations of physics are at the same time overdetermined and 
underdetermined, due to gauge symmetries, see [ 17 ], [ 18 ], [ 19 ], [ 20 ]. The variational origin 
of equations also lead us to assume that equations live in the dual space of field deformations, 
which leads us to a framework in which the metric on fields is not needed to define symmetric 
hyperbolic systems. 

Further investigations are needed to generalise this to more general models of interest for (fun¬ 
damental) physics. It seems possible that integrability is generically equivalent to Hamiltonian 
formulation, though in a completely Lagrangian setting. 

We also found a more general class of coordinate systems than harmonic coordinates for 
which the system becomes symmetric. This is not that important in vacuum gravity since 
harmonic coordinates always exist and they are sufficient to solve the CP. However, when 
gravity is coupled to matter fields, matter field equations may need to be symmetrised as well. 
Unfortunately, the choice of spatial coordinates in a game which can be played only once since 
matter equations are coupled to gravity equations. One would need a coordinate system in 
which the whole system is symmetric, while the symmetrisation of matter equations depends 
on the matter-gravity coupling. Having more coordinate systems in which gravity becomes 
symmetric hyperbolic may help when coupled with some matter fields. We still do not have 
examples of matter fields which can be solved in this way, though it is clear that harmonic 
coordinates play no distinguished role within the class of coordinates which turn the system in 
a symmetric hyperbolic system. 
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Appendix A. Projectors 

Let us define the projectors for the ADM decomposition. Due to the immersion of the spacelike 
hypersurface S (t: E ^ M), we have a normal covector: 

u = ^e^^^dAx'"dBx‘'dcx^€f,i,padx'" (A.l) 

and, due to the metric structure, its normal vector: 


u = g^''upd„. (A.2) 

Since ft is not lightlike, we have its normal unit vector n, obtained by normalization. 

With n we can define a basis adapted to the foliation: let us take x = i(k) € M we can take 
a quadruple of vectors: (n, ca) where ba is a basis in T^E. 

This set has the propriety that each ca is orthogonal to n, so it is a basis for . Also we 

have that the norm of n is —1, i.e.: 

g{n,n) = -l. (A.3) 

Now, we can define some maps which allow us to decompone each geometrical object (like 
tensors, metrics and so on). These are defined as follows: 

(j'i = + n^ni,. (A.4) 

and it is easy to see that they are idempotents, so they are projectors. Namely, we have: 

<< = (^-5) 

Once projectors are defined we can decompose vectors and covectors (and, in general, tensor) 
in tangent and normal part. 

Theorem: For all v = v^dp € T^(k)^ there exists a decomposition: 


V = Vjl + v± 


with U|| G TfcE. This decomposition is unique. 

Proof: 


V = v>^dp = v^dp = v^5';d,. 


We can obtain (by using (A.4)): 


(A. 6) 


so that: 


Sit = att 


— n^riv 


V = v^(a'( - rfripOdv = - v^rAupdi, = v^a'td^ + {-v^npOrA 


and we have: 


V = V|| + v± 


(A.7) 
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wit W|| = v^a’i^di, and wx = —It is easy to see that g{v±,eA) = gin,v\\) =0. ■ 

A similar result holds for the covectors. Let a a covector on (i(S)), namely: 

a = a,dx^ e (A.8) 

Then there exist a decomposition of a (parallel and perpendicular part) and a covector [3 on 
S such that: 


6*(a||)=/3eT*(E). (A.9) 

Let now g = g^^(x)dx^ 0 dx’' be a metric on i(E) and we can write its tangent part: 

9\\ = gap{x)<7la'^dx^dx''. (A.IO) 


Now, we can take the pullback along the immersion: 

L*{g\\) =gap{x{k))<j'^a^dAX^dBx''dk^ ® dk^ = 
gapiS'^ + + n'^ni,)dAX^dBx''dk^ (g) dk^ = 

{g^v + n^inv)dAX^dBx''dk^ ® dk^ = g^i^dAX^dBx''dk^ ® dk^ = 
yAsdk^ ( 8 > dk^ = 7 


(All) 


This means that the parallel part of the metric restricts to the induced metric over S. 


Appendix B. 

Let us here consider an example for the situation discussed in Section 2. Fix two fields and 
and one spatial coordinate on the space E (i.e. dim(M) =2). In this simple case equation 
(2.13) reads as 


fffil 

/3i2 

0 

0 

0 



(y^\ 


/o 

0 

0 

0 

0 





1^22 

0 

0 

0 

0 




0 

0 

0 

0 

0 

0 


y^ 

0 

0 

0 

0 

ail 

<221 

ai2 

0:22 

0 

0 

0 

0 

5o 


- 

0 

0 

0 

0 

a\i 

ah 

0^12 

0^22 

aj} 

ah 

a}? 

di 


0 

0 

0 

0 

PW 

p\l 




0 

0 

Pll 

Pll 

0 

0 



\ 0 

0 

0 

0 


Pl\) 


\vjj 



0 


Pll 

0 

0 y 


\vV 


(B.l) 


from which one sees that in order to have a symmetric hyperbolic system one needs to verify 
condition (2.14). 

Let us stress that the symmetry of the system relies on a suitable ordering of fields and 
equations. This is acceptable since symmetric hyperbolic form is sufficient, not a necessary, 
condition for having existence and uniqueness of solutions. 
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